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Abstract 



We extrapolate lattice calculations of electric charge radii of the spin- 1/2 
baryon octet to the physical regime. The extrapolation procedure incor- 
porates chiral perturbation theory and heavy quark effective theory in the 
appropriate limits. In particular, this procedure includes the non-analytic, 
logarithmic terms from pion loops. The electric charge radii of the nucleons 
and S~ obtained from the chiral extrapolations agree "well "with experimental 
data. We make predictions for the charge radii of the remaining baryons in 
anticipation of future experimental measurements. 
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I. INTRODUCTION 



Lattice QCD is so far the most successful method of calculating hadronic observables from 
the theory of QCD. However, computational limitations mean that hadronic observables are 
calculated on the lattice at quark masses larger than their physical values. Hence results 
from lattice simulations cannot be directly compared with experimental data. Although, 
with improvements in actions, algorithms and computer speed, future lattice calculations 
will be performed much closer to the physical regime, these improvements will proceed over 
many years. Therefore, to make sense of any lattice results produced to date, and to compare 
them with experiment, one must understand how to extrapolate lattice results, obtained at 
large quark masses, to the physical world. 

One of the difficulties with calculating hadronic observables at heavy quark masses on 
the lattice is that chiral perturbation theory cannot be applied in this quark mass regime. 
However, chiral expansions of hadronic observables contain important non-analytic terms 
as a function of the quark mass, rriq (or equivalently of m^, as nig oc in this range). 
It is vital that this non-analytic behaviour is included in any extrapolation to the physical 
regime §. 

The chiral expansion of the squared electric charge radius of a spin-1/2 octet baryon 
includes non-analytic behaviour in the form of logarithmic terms in (or nig). To extrap- 
olate the lattice results for electric charge radii we incorporate these logarithmic terms in 
our extrapolation formulae, while ensuring that the correct heavy quark behaviour is also 
maintained. A similar approach has been successful in explaining why lattice simulations 
of pion and proton charge radii are similar in size, while experimental measurements reveal 
a significant difference. The dramatic differences in the chiral behaviour of the pion and 
proton charge radii account for the similarity of the lattice results at moderately heavy pion 
masses, while allowing good agreement with experiment at the physical pion mass. This il- 
lustrates the importance of including meson cloud effects in extrapolations of lattice results 
to the physical regime. In this paper, we improve the formalism of chiral extrapolations by 
incorporating both chiral symmetry and heavy quark effective theory. 

Our lattice "data" for the electric charge radii is taken from the calculations of Ref . . 
The contributions to the charge radii from individual quark flavours are also given there. 
These results are the only available lattice calculations of the electric charge radii of spin-1/2 
octet baryons. For each baryon, the electric charge radius was calculated at three different 
quark masses, corresponding to rather heavy pion masses (all exceeding 600 MeV). Here we 
extrapolate these results as functions of the squared pion mass, m^, to obtain predictions 
for the charge radii at the physical pion mass (139 MeV). Because the lattice calculations 
are quenched, we expect that there are errors in the lattice data which we have been unable 
to take into account. However, we expect that errors from the quenching approximation 
will be rather small. At the quark masses considered on the lattice, the dominant effect is a 
simple renormalization of the strong coupling constant, accounted for in setting the lattice 
spacing scale. When results from unquenched simulations become available, the formalism 
presented here may be readily applied. 
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II. EXTRAPOLATIONS 



Dynamical breaking of chiral symmetry in the QCD lagrangian results in the formation 
of an octet of (pseudo-) Goldstone bosons. Goldstone boson loops give rise to significant 
non-analytic behaviour in hadronic observables, such as (r-^) and magnetic moments, as a 
function of the quark mass, rrig. Using an expansion about the chiral SU(3) limit gives the 
following expression for the squared electric charge radius, (r-^), of a spin-1/2 octet baryon 
(labelled by i) 



X=TT,K 



Here is the pion decay constant (93 MeV) and A is the scale of the dimensional regular- 
ization. (The value of 7, is clearly correlated with the choice of A.) Unlike SU(2)-flavour 
symmetry, SU(3) is significantly broken in the physical world, with the strange quark mass 
the same order of magnitude as Aqcd for low-energy phenomenology. The squared kaon 
mass exceeds the squared pion mass by over an order of magnitude. Given that the source 
of the meson cloud associated with a baryon is of a finite size, one might anticipate that the 
role of the kaon cloud will be suppressed away from the SU(3) chiral limit The form 
factor describing the finite size of the kaon source will act to suppress kaon loop effects like 
nix'^ at large mx (comparable to mx)- This was demonstrated within a particular chiral 
quark model for the nucleon magnetic moments in Ref. [Q. Despite the model-dependence 
associated with the form factors, the lattice results themselves are very slowly varying func- 
tions of mx at values of the order rriK or higher, thus supporting the general conclusion. 
Therefore kaon loop effects are expected to be small and slowly varying as a function of 
rriq. Hence we do not explicitly include kaon contributions in our extrapolation formulae. 
Conversely, since varies rapidly with rUg, the leading non- analytic behaviour in must 
be included explicitly in an extrapolation to the chiral regime. 

To extrapolate the lattice calculations of the electric charge radii of the spin-1/2 octet 
baryons, we consider two distinct fitting procedures. Both these extrapolation schemes 
satisfy the constraints of chiral perturbation theory and heavy quark effective theory. The 
first extrapolation procedure we investigate is given simply by the formula 

2\ Ci + Xiln(m^/A) 



/ 1 + C2m^ 

where (r,^) are the lattice QCD results (at several values of m^) extracted from Ref. 0, 
ci and C2 are fit parameters chosen to best fit these results, Xi (corresponding to the i^^ 
baryon) is fixed (model-independently) by chiral perturbation theory and A, which is directly 
correlated with ci, is fixed at 1 GeV. This extrapolation procedure is not feasible for the 
neutral baryons because (r^^) ^ as becomes large and thus sensitivity to the C2 fit 
parameter is lost. In order to extrapolate the neutral baryon charge radii results we consider 
a second extrapolation procedure, focusing on individual quark "sector" (or quark flavour) 
contributions, as discussed below. 
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Clearly the extrapolation formula given in Eq. (|]) builds in the correct chiral behaviour 
in the SU(2) limit, since in the limit m^r — it can be expanded as follows 



Ti ) = Ci + Xi In(m^) - C2m^ + . . . (3) 

(Recall that the scale A, in Eq. (2), has been set to 1 GeV. This choice is also implicit in 
Eqs. (6)-(10), below, where ttItt in the logarithm must be in GeV.) This agrees with the 
chiral SU(2) expansion of the squared electric charge radius (see Eq. (||)), provided we fix 
the coefficient Xi to Qaf^^ / {A7rf„)'^ for the i^^ baryon. The coefficients af^^ and Xi ^-^e given 
in Table |. 

In the large rriq limit, we expect that the quarks behave non-relativistically, and the 
squared charge radius falls off as as it does in non-relativistic quantum mechanics. In 
the region where rriq is very large, oc m^, and hence we require that 

r^') oc ^, (4) 

as becomes extremely large. This is clearly satisfied by our first extrapolation formula, 
Eq. (^, since the logarithm is very slowly varying. 

In the second extrapolation procedure we deal separately with the individual quark 
sector contributions to the baryon charge radii. For example, in the case of the nucleons, we 
extrapolate the up and down sector contributions separately. For the hyperons the strange 
and light sector results are extrapolated separately. This avoids the problem encountered 
with the neutral baryons which was mentioned previously, because now all the quantities 
being extrapolated are charged, even if the overall charge on the baryon is zero. This 
separation is valid because the squared electric charge radius can be decomposed as 

r^)= E e,(rl^^^), (5) 

q=u,d,s 



where (rl y is the contribution from the q quark sector and Cg is the charge of this 
quark sector. Therefore, provided that the extrapolation formulae from each sector add so 
that the chiral and heavy quark limits of the sum are in agreement with Eqs. (|I]) and (^ 
respectively, this method contains the same physics as the first method, but simply makes 
use of the extra information contained in the individual quark sector results. Not only does 
this second extrapolation procedure solve the neutral baryon extrapolation difficulty, it also 
provides predictions for individual quark sector radii, which will be resolved at Jefferson Lab 
for the nucleon |]12|, and perhaps future accelerator facilities for hyperons. 



Isolation of the individual quark sector contributions to the charge radii is relatively 
straightforward from the theoretical point of view. For example, to isolate the w-sector 
contribution to the charge radius of the proton one simply sets the cZ-quark charge to zero 
and calculates the proton charge radius as if only the u quark carried charge. In the chiral 
expansion of the proton charge radius, the coefficient of the logarithm, x^\ originates from 
the pion loop, p — > n+7r+ (see Fig. and includes the charge of the pion cloud. Therefore, to 
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FIG. 1. Schematic illustration of the pion loop which produces the leading non-analytic con- 
tribution to the proton charge radius. 



extrapolate the w-sector contribution to the proton charge radius, the appropriate coefficient 
of the logarithm is ^X^\ since the pion now carries charge +2/3. Thus we extrapolate the 
M-sector results of the proton according to 

2e„(.(-)^) = £llMii!^. (6) 



where x^'^ is the full chiral coefficient of the proton, given in Table | and (r^^'^J is the 
squared charge radius of a single w-quark of unit charge. Similarly, the ci-sector results are 
extrapolated according to 



eJrif)^\ = ^±^^^EppLi, (7) 



.(d) 2\ _ c/+ 3Xpln(m ^ 
P I 1 + C2'ml 

where the factor of 1/3 originates from the d contribution to the pion cloud. Clearly adding 
the left hand sides of Eqs. (|^) and yields the full expression for (j'f,'^- In the chiral 
limit the right hand sides add so that the correct chiral form for (j'p'^, given in Eq. (j^), is 
retained. The sum of Eqs. (|^) and (J^) also obeys the correct heavy quark behaviour, given in 
Eq. (^. Since the parameters C2 and C2' are not necessarily the same, the individual quark 
sector extrapolation formulae (Eqs. (§) and (0)) cannot be added directly to give Eq. @. 
Therefore, in general we do not expect the two extrapolation procedures to give exactly the 
same results. For the charged baryons this may be used to help quantify the systematic 
error of the approach. For the neutron, the analogous extrapolation functions are given by 

Ju)2\ _ Cl + fxl"^ H^T.) , . 



for the M-quark sector, and 



2e.(r(^)^) = £L4M^_i!H!^, (9) 



1 + C2 'ml 

for the d-quark sector results, where Xn"^ is given in Table |. 
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We now consider extrapolating the hyperon charge radii results using the second extrap- 
olation procedure, i.e., extrapolating the strange and light quark sector results separately. In 
extrapolating the strange quark sector the charges of the light quarks are set to zero. Since 
the logarithmic term considered here originates from pionic corrections to the charge radius 
(and pions do not contain strange quarks), the coefficient of the logarithm in the strange 
sector extrapolation will be zero. Similarly, for the light sector extrapolation, strange quarks 
do not carry charge and hence the coefficient of the logarithm in this extrapolation will be 
the full coefficient, x!'^''- This results in the following extrapolation formulae for the hyperon 
quark sector contributions 

^/ ,0.V c. + x.lu(m.) 

\ / 1 + C2m^ 

and 



W2\ _ Ci 



1 + C2 'ml 



(11) 



where i runs over the hyperons only, and I corresponds to the light-quark {u and/or d) 
sector. Since the strange quark mass is held fixed in the light quark mass extrapolation, any 
variation in the strange quark sector is purely an environment effect from the surrounding 
light quarks. As such, the functional form for the strange quark sector is constrained by 
neither leading order chiral perturbation theory nor heavy quark effective theory. As we 
shall see, C2' is small and negative for each hyperon, which suggests that a simple linear 
ansatz for the strange quark sector extrapolation could also have been used. 



III. RESULTS 

The lattice QCD simulations were performed on a 24 x 12 x 12 x 24 periodic lattice 
using standard Wilson actions at /5 = 5.9. Dirichlet boundary conditions were used for 
fermions in the time direction. Twenty-eight quenched gauge configurations were generated 
by the Cabibbo-Marinari pseudo-heat-bath method. The conserved vector current was 
derived from the Wilson fermion action via the Noether procedure. The associated lattice 
Ward identity protects this vector current from renormalization. The radii were produced 
by fitting the electric form factor to dipole and monopole forms, allowing a charge radius 
to be extracted in each case. Since it is known from experiment that the dipole form is 
more suitable for parameterizing the electric form factor, we consider only the dipole results 
here. Statistical uncertainties in the lattice simulation results are calculated in a third-order, 
single elimination jackknife ]T^ , pD| . Further details may be found in Ref. [^]. 



The extrapolations of lattice calculations for the charge radii of the spin-1/2 baryon octet 
are shown in Figs. Extrapolations of baryon charge radii results, performed according 
to Eq. (0), are indicated by the solid lines, where the full circles (•) represent the baryon 
charge radii from lattice QCD and the extrapolated value at = 139 MeV. The individual 
quark sector extrapolations are shown by the dashed and dot-dashed lines, and the baryon 
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TABLE I. Baryon electric charge radii and the quark sector contributions. The latter are 
defined on the left hand sides of Eqs. (|^-(11). One- loop corrected estimates of af^ (in Eq. (jl|)) 
and Xi (ill units of fm^) for each octet baryon are indicated. For each extrapolation, the fit 
parameters, ci and C2, and the predicted value of (r^) at the physical pion mass are reported. 
Asterisks denote the squared charge radii reconstructed from the sum of separate quark sector 
extrapolations. (The units are such that the pion mass is in GeV and the squared charge radius 
in fm^.) 7 
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FIG. 2. Fits to lattice results for the squared electric charge radius of the proton. Fits to the 
individual quark sector results are also shown. The u-quark sector results are indicated by open 
triangles and the d-quark sector results by open squares. Physical values predicted by the fits 
are indicated at the physical pion mass, where the full circle denotes the result predicted from 
the first extrapolation procedure and the full square denotes the baryon radius reconstructed from 
the quark sector extrapolations (see text). (N.B. The latter values are actually so close as to be 
indistinguishable on the graph.) The experimental value is denoted by an asterisk. 



charge radius predicted by this method is indicated by a full square (■). Experimental 
measurements are indicated at the physical pion mass by an asterisk (★). Note that for 
the charged baryons, two extrapolation schemes and two corresponding predicted physical 
values are shown, whereas (for the reasons explained in section ^ only one extrapolation 
procedure is shown for each neutral baryon. 



In the case of the proton, the two extrapolated results agree very well with the exper- 
imental measurement. It can be seen that a traditional linear extrapolation in would 
significantly underestimate the experimental result. Similarly the predicted charge radius 
for the neutron (produced by separate extrapolations of the u- and rf-sector results) agrees 
with the experimental data significantly better than the prediction from a conventional lin- 
ear extrapolation in m^. Both predicted values for E~ also agree very well with the two 
experimental measurements. These baryons are currently the only baryons of the spin-1/2 
octet whose electric charge radius has been measured. 

From simple quark model arguments where the heavier strange quark has a smaller 
distribution than the light quarks, we expect the hierarchy of electric charge radii of charged 
octet baryons to be given as follows 



> 



> 



> 



(12) 



Clearly the results of the lattice extrapolations shown in Table | are in qualitative agreement 
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FIG. 3. Fits to lattice results for the quark sector contributions to the squared electric charge 
radius of the neutron. All symbols are defined in the caption below Fig. 0. 



with this expectation. Indeed, in the regime of the actual lattice data (m^r > 600 MeV) 
the argument is even quantitative. However, as the chiral limit (and physical pion mass) 
is approached, the simple quark model description is no longer adequate and chiral physics 
gives rise to dramatic effects. For example, in the extrapolation of the (i-quark sector of 
the proton (Fig. chiral effects mean that the ci-quark sector can actually make a positive 
contribution to the charge radius, via the d contribution in 7r+ - even though the total 
contribution is negative at the physical pion mass. This behaviour is not anticipated by the 
simple quark model. 

For the neutral baryons the sign of the squared charge radius is important. In the 
neutron, the two d quarks are most likely to be found in a spin 1 configuration, where they 
will undergo hyperfine repulsion. This leads to a small, negative charge radius squared. 
However, as one approaches the chiral limit, spontaneous chiral symmetry breaking, in 
particular the process n Ptt^, which carries d-quarks to larger radii and screens the 
M-quark contribution via the u in the vr", leads to an enhancement of the negative charge 
radius. The remaining neutral baryons, A and have a positive squared charge radius. 
This is because in each case the strange quark distribution is more localized than the light 
quark charge distribution (due to the larger mass of the strange quark). Therefore on average 
the light quark charge distribution occurs at a larger radius, resulting in a positive charge 
radius (since the light quark charge is positive in each case). 

The lattice results used here were calculated with a strange quark mass of approximately 
250 MeV 1^. This is much larger than the physical strange quark mass of 115 ± 8 MeV 
at a scale 2 GeV, taken from a careful analysis of QCD sum rules for r decay [|l^. In an 
earlier study of lattice results for octet baryon magnetic moments ||^ (where the results 
were extracted from the same lattice simulation 0), it was found that the heavy strange 
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FIG. 4. Fits to lattice results of the squared electric charge radius of the S"*". Fits to the 
individual quark sector results are also shown. The strange quark sector results are indicated by 
open diamonds. All other symbols are defined in the caption below Fig. 0. 
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FIG. 5. Fits to lattice results of the squared electric charge radius of the YT . Fits to the 
individual quark sector results are also shown. All symbols are defined in the captions below 
Figs. I and |. 
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quark had a significant effect on the predictions of the S moments. Here we expect that 
the heavy strange quark should also have some effect on the S charge radii. With a strange 
quark mass closer to the physical mass the strange quark contribution would be increased. 
This would result in a lower predicted charge radius for the S° and a larger (in magnitude) 
charge radius for the H~. In the absence of experimental measurements we will not attempt 
to correct for the effect of the strange quark mass here. 

As we see from Table |, the extrapolated mean square charge radii obtained from both 
extrapolation procedures agree quite well for each charged octet baryon. For the proton, 
and 5~ the reconstructed values completely cover the result from the original extrapolations 
of Eq. (0). In the case of the E"*" the two values overlap only on the error bars. This is due 
to the small variation in the strange quark contribution (which is caused by an environment 
effect). When this environment effect is included in the baryon charge radius, the magnitude 
of the slope is increased, resulting in a larger charge radius after extrapolation. 

In turning the dimensionless masses calculated on the lattice to physical units, the lattice 
spacing, a, was set in the traditional manner by fixing the nucleon mass, obtained by a naive 
linear extrapolation in rri^, equal to the observed mass. Of course, such a linear extrapolation 
is known to be inconsistent with chiral symmetry. We have checked that applying a more 
consistent chiral extrapolation would systematically lower values of (r^) obtained for the 
charged octet by of the order 15%. (The effect on neutral baryons is much smaller.) On 
the other hand, the data which we are forced to use is quenched data which omits some 
pion corrections. Although these are expected to be suppressed at the large values of 
for which the data is available (c.f. Ref. [^]), the associated systematic error would tend to 
increase the calculated values of (r^), perhaps by 5-10%. Rather than attempt to repair these 
deficiencies in the present data, we feel it would be more reasonable to simply accept that 
there is an additional systematic error of the order 15% associated with the extrapolated 
values shown in Table |. Adding this systematic in quadrature means that the values in 
Table | would become, for example, (r?^^ = 0.68 ± 0.14 fm^, {r'^^_) = -0.54 ± 0.09 fm^, 
(r?2+) = 0.77 ± 0.14 fm^. We look forward to repeating our analysis with unquenched data 
at lower quark mass, which is the best way to overcome these problems. 

IV. CONCLUSION 

In this paper we have investigated two methods of extrapolating lattice results for the 
electric charge radii of octet baryons to the physical regime. These procedures build in the 
correct leading non-analytic behaviour of the electric charge radii in the chiral limit, as well 
as the correct heavy quark behaviour. Both extrapolation procedures were performed for 
the charged octet baryons, with the predicted values agreeing very well. The extrapolation 
formulae seem to be very successful, as good agreement with experiment was obtained for 
the nucleons and the S~. We await further experimental measurements of the baryon charge 
radii in order to test our predictions. In the future we hope to perform similar extrapolations 
of electric charge radius lattice results calculated with a more realistic strange quark mass 
and eventually with lighter, dynamical u and d quarks. 
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